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Abstract: By means of fluctuational electrodynamics, we calculate 

radiative heat flux between two planar materials respectively made of SiC 
and Si02. More specifically, we focus on a first (direct) situation where one 
of the two materials (for example SiC) is at ambient temperature whereas 
the second material is at a higher one, then we study a second (reverse) 
situation where the material temperatures are inverted. When the two fluxes 
corresponding to the two situations are different, the materials are said 
to exhibit thermal rectification, a property with potential applications in 
thermal regulation. Rectification variations with temperature and separation 
distance are reported here. Calculations are performed using material 
optical data experimentally determined by Fourier transform emission 
spectrometry of heated materials between ambient temperature (around 300 
K) and 1480 K. It is shown that rectification is much more important in 
the near-field domain, i.e. at separation distances smaller than the thermal 
wavelength. In addition, we see that the larger is the temperature difference, 
the larger is rectification. Large rectification is finally interpreted due to 
a weakening of the SiC surface polariton when temperature increases, a 
weakening which affects much less Si 02 resonances. 
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1. Introduction 

Since the seminal works of Rytov 12 and Polder and Van Hove ||2l, it has been known that 
radiative heat transfer (RHT) can be much larger than classical Planck’s prediction when the 
separation distance between two solid bodies are small compared to the thermal wavelength at 
which heat radiation is exchanged. Enhancement is particularly large in the case of dielectrics 
supporting phonon-polariton such as SiC or Si02 IS that have probably been the most stud¬ 
ied materials in the framework of near-held heat transfer. In the case of materials supporting 
polaritons, these non radiative modes can couple so that heat is transferred in vacuum through 
photon tunneling. However, the magnitude of the enhancement strongly depends on the ampli¬ 
tude and the width of the resonance in the material optical response which can be affected by its 
temperature HHJl. Indeed, when the temperature increases, one knows that phonons in general 
including phonon-polaritons are affected by several processes like phonon-phonon interactions 
so that theses modes see their lifetime decreases ||6l. But phonons in different materials are not 
equally affected by temperature changes. 

A question which naturally arises is whether the heat flux between a material at ambient tem¬ 
perature and another one different at a higher temperature will be the same when temperatures 
are inverted. If it is not the case, the two materials couple is said to exhibit thermal rectihcation, 
a thermal analog of electrical rectihcation better known as the diode. 

Thermal rectihcation has recently received a growing interest due to the emergence of ther¬ 
mal How management needs related to limited energy and global warming issues. If thermal 
equivalent of the diode and the transistor would be designed and fabricated, this could open the 
way to new passive temperature regulation and thermal circuits with no needs of electronics 
leading to very reliable thermal management. 

In the last decade, thermal conductors exhibiting thermal rectihcation have been proposed. In 
these devices, heat carriers hux is made asymmetric by nanostructuring the material or by guid¬ 
ing the carriers Hux differently whether they are moving in one direction or the other ||7l4T5]l . 
These physical ideas can be used to conceive and fabricate a thermal transistor leading to the 
possibility of logical thermal circuits Gsiini. More recently, radiative thermal rectihers have 
been proposed ifTStESll as well as radiative thermal transistors II23II24II . Very promising devices 
have been imagined based on materials exhibiting phase change transition such as Vanadium 
dioxyde (VO 2 ) ll^l^l25]l . superconductors f22\ or thermochromic materials ll26ll . The goal 
of this article is to show that radiative thermal rectihcation can also be achieved with the most 
studied polaritonic supported materials like SiC and Si02. 

We hrst present optical properties measurements showing how both SiC and Si02 dielectric 
functions change with temperature. These measurements performed using Fourier Transform 
Emission Spectroscopy allow us to make near-held radiative heat transfer calculations taking 
into account the actual material dependance with temperature. Similar near-held calculation 
taking into account the temperature dependance of optical properties have already been taken 
in the past lIZTl with measurement on 3-C SiC of ll28ll but in a different situation between SiC 
and metallic coated-SiC. Here, we calculate heat transfer between two planar interfaces of SiC 
and Si02. We show that these two materials exhibit almost no thermal rectihcation when the 
exchange occurs in the far-held i.e. when the main contribution to the hux comes from the 
propagative waves. On the contrary, at distances below 100 nm, we show that rectihcation ratio 
drastically increases in order to reach values as high as the highest predicted and measured 
ones in the literature. We hnally show that rectihcation is mainly due to the SiC resonance 
attenuation at high temperature. 




2. SiC and Si02 optical datas 

In this section, we report SiC and Si02 optical properties temperature dependance with tem¬ 
perature. The dielectric functions of SiC and vitreous Si02 have been extracted from measure¬ 
ments of their spectral emittance at near normal incidence. The SiC sample provided by MTI 
Corporation is a single crystal with a (0001) orientation and two epi polished sides. The stacking 
sequence is of 6H type and its electric resistivity lies between 0.02 and 0.2 ohm.cm. The Si02 
material is a polished plate of vitreous silica with low OH content (< 2Qppm). The emittance 
spectra were acquired with a homemade spectrometer built around two Bruker infrared spec¬ 
trometers, a Vertex 80v working under vacuum and a Vertex 70 purged with dry air (see rtSOl l. 
The first and second devices have been used to probe the far- and mid-infrared ranges respec¬ 
tively, at a spectral resolution of 4 cm^*. In order to avoid parasite fluxes during the measure¬ 
ments, a CO 2 laser has been used to heat the samples up to 1500 K. The retrieval of the optical 
functions has been performed in a second step by fitting the experimental data with adequate 
physical expressions. A semi-quantum dielectric function model, well adapted for crystalline 
media, was selected to reproduce the infrared response of the silicon carbide sample ll^ . The 
disordered nature of vitreous silica needed instead the use of a causal Gaussian dielectric func¬ 
tion model able to take account for the inhomogeneous broadening of the absorption bands of 
glasses ED. The whole process allowed to obtain the temperature dependences of the opti¬ 
cal properties for both materials in a wide range, the results are reported in Fig. [D For Si02, 
we note that the resonance around 2x10*^ rad s^* slightly shifts and broadens when the tem¬ 
perature increases. For SiC, we note a similar behavior (shift and broadening) but even more 
pronounced for the resonance around 1.5 x lO'^ rad s^*. This is due to the fact that anharmonic- 
ity effects responsible of the phonon-phonon interaction increases with the temperature ja.we 
will see that this behavior is responsible for the important rectification between Si02 and SiC 
in the near-held. Note however that the second resonance in Si02 does not change a lot with 
temperature as long as temperature is below the vitreous transition which occur around 1500 
K E21- This is probably related to the very weak Si02 thermal dilatation. 

3. Radiative heat transfer between SiC and Si02 

We now come to the calculation of the RHT between two planar interfaces made of SiC 
(medium 1) and Si02 (medium2). RHT between two planar interfaces can easily be calcu¬ 
lated using Huctuational electrodynamics formalism. Heat Hux appears as a semi-analytical 
expressions which is nothing else than the summation of individual plane wave contributions to 
the heat transfer. These plane waves are labelled by a triplet {(0,K,a) where co is the angular 
frequency of the wave, K is the wavector component parallel to the interfaces whereas a is 
the wave polarization (s or p). Let us call Ti temperature of medium 1 and T 2 temperature of 
medium 2. Heat flux expression between media 1 and 2 reads iBlI^ 



( 1 ) 


where 0(0,7’) = h(o/[exp[h(o/khT] — 1] is the mean energy of a photon at temperature T. 
In the preceding expression, t‘^{co,K) appears as a transmission coefficient for a plane wave 
{co,K,a) between medium 1 and medium 2 or medium 2 and medium 1. Note that expres¬ 
sions of transmission coefficient are different whether the wave is propagative {K < (o/c or 
evanescent {K > co/c). If we call medium 3, the medium separating the two other media and 
introducing Fresnel reflection coefficients, the expression for the transmission coefficient reads 
for propagating waves 
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Fig. 1. Measured real and Imaginary part of the dielectric function e variation with angular 
frequency co for different temperatures. Top 2 figures ; Si02. Bottom 2 figures : SiC. 
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In these expressions, 3 denotes the imaginary part, 7 ; is the perpendicular wave vector compo¬ 
nent in medium i which for non magnetic materials reads Jj = The usual fresnel 

reflection coefficient are given by r\j = {ji — Yj)/{ji + Yj) and r?- = {EjYi — eiYj)/{£jYi + ^iYj) 
when one considers local materials which is the case at nanometric separation distance II34II35I . 
The flux expression only depends on the material temperatures, their local optical response and 
the separation distance between the bodies. Note that the symmetry of Eq. ([T]i implies that recti¬ 
fication will occur if and only if the two materials are different and have their optical properties 
that change with temperature in a different manner. 

Optical properties variations with temperature have rarely been taken into account, mostly 
in the case of phase change materials. Here, our main goal is to perform these calculations in 
the case of two very well studied materials but for which optical properties variations with tem¬ 
perature are taken into account. The input properties in our calculation are the optical properties 
presented in the preceding section. It can be noted in these measurements that they have been 
taken at close temperatures for both materials but not exactly at the same one. As explained 
in the introduction, our goal is also to explore how a system constituted of SiC and Si 02 can 
exhibit radiative thermal rectification. Thus, a difficulty arises since rectification compares a 
situation where the materials are at two different temperatures with a second situation where 
the temperatures considered before have been inverted. For example, measurements at ambient 
temperature for SiC have been performed at 295 K whereas they have been performed at 298 K 
for Si 02 . In our calculation, we choosed ambient temperature to be 297 K and we assumed that 
optical properties at this temperature were the ones at 295 K for SiC and 298K for Si 02 . Other 
temperatures considered in the paper are 471 K, 671 K, 981 K, 1102 K and 1470 K. To perform 
heat flux calculations, optical properties for SiC and Si 02 have been respectively taken at 464 
K and 479 K, 686 K and 672 K, 996 K and 966 K, 1105 K and 1098 K, 1460 K and 1480 K 
(See Fig.[T])- 

We report in Fig. |2] heat flux variations as a function of the separation distance between the 
two materials. We compare the situation where Si 02 is at 297 K and SiC at 1470 K with the 
situation where the temperatures are inverted. We note that for separation distances larger than 
a few microns the two fluxes are constant and almost superimposed each other. When Si 02 is at 
ambient temperature, the value of the heat flux is 1.872x 10^ W m^^ whereas it is 1.896x 10^ W 
m^^ in the reverse case. These two values are around 71% of the heat exchanged between two 
blackbodies for which Fresnel reflection factors are equal to 0 in the fluxes expression. When 
the separation distance is decreased, the heat flux increases for both situations. For distances 
larger than 200 nm, heat fluxes are very close in both situations. For distances smaller than 200 
nm, we see that the situation where SiC remains at ambient temperature exhibits a higher flux 
than the one corresponding to the reverse one. 

We define rectification R as the relative variation of the heat flux 1120112111 in the two situations 
so that 

^ ^ 1^102(^1,72) - ^102(72,7^1)! 

Max[^io2(7’i, 72), 9102(72, T’l)] 

With such a definition, R = 0 means the the material exhibit no rectification whereas R = 1 
correspond to the case of a perfect thermal diode. In Fig. [3 we report the variation of the 
rectification as a function of the separation distance between the two planar interfaces made of 
SiC and Si 02 for different temperatures. One material is always at 297 K whereas the second 
one is taken at 471 K, 671 K, 981 K, 1102 K or 1470 K. We note a very similar behavior for 
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Fig. 2. Computed radiative heat transfer between a plane interface of SiC and a second one 
constituted of Si02 versus their separation distances. Two situations are compared. In the 
first one (plain) Si02 is at 297 K and SiC at 1470 K whereas in the second one (dashed) 
SiC is at 297 K and Si02 is at 1470 K. 



Fig. 3. Computed rectification variations as a function of the separation distance between 
two planar interfaces made of SiC and Si02 when one material is at 297 K and the second 
one at 471 K, 671 K, 981 K, 1102 K or 1470 K. 
















all the cases presented. For distances larger than 200 nm rectification remains very low. It goes 
from 10^^ for the smallest temperature difference to 10^^ for the largest. One can therefore 
say that SiC and Si02 exhibit very low radiative thermal rectification in the far field and more 
generally for separation distances larger than 200 nm. Below 200 nm, rectification increases 
drastically to reach values up to 0.7 for a temperature of 1470 K. Rectification saturates when 
the separation distance goes below 10 nm. Note also that the higher temperature difference is, 
the higher is the rectification in the near-field. 

4. Discussion 

Let us now shed some light on the physical reason that is responsible of the behavior of the 
rectification between SiC and Si02 with both distance and temperature. To do so, we plot 
in Fig. |4] the variation of the spectral radiative heat transfer flux as a function of the angular 
frequency for the case where one material is at 297 K and the other one at 1470 K. The spectral 
flux is represented for 4 different separation distances 1 jim, 100 nm, 10 nm and Inm. At 
large distance (here 1 /rm), the spectral flux is very similar to the one representative of the 
exchange between two black bodies except that there exist some frequencies where one can 
identify dips. These dips correspond to the frequencies where Si02 and SiC exhibit resonances. 
Indeed, Si02 exhibit two resonances and SiC one resonance that had been identified in the 
first section. At these frequencies Si02 and SiC are very reflective so that transfer is greatly 
reduced. However, the range of frequencies where materials are low absorbent is much smaller 
than the one where both materials almost behave as blackbodies. This explains why at large 
distance when this spectral flux is integrated over all frequencies, heat flux reaches 70% of the 
one between two blackbodies. If we examine carefully the region where the spectral flux is 
much lower than between blackbodies, one see that these are the regions where the spectral 
fluxes differ in the direct and the reverse situations. We notice that the depth of the dips is less 
pronounced when the material concerned is at high temperature. For example dips for Si02 
(around 9xl0'^ rad s^* and 2x 10*^ rad s^') are less pronounced in the plain curve where 
Si02 is at high temperature. In the same way, the dip corresponding to the SiC resonance 
around 1.5 x 10*^ rad s^' is smaller in the dashed curve corresponding to the case where SiC 
is at high temperature. Despite of these differences, as most of the contribution comes from the 
region where the spectral fluxes are equal, rectification in the far-field is very small and even 
smaller at lower temperature when the dielectric function is still close to the one at ambient 
temperature. 

When the distance is reduced, the situation drastically changes. The main contributions to 
the spectral fluxes mainly comes from the frequencies where the resonances are present. This 
is a well known behavior when one is treating materials supporting phonon polaritons such as 
SiC and Si02. Indeed, the heat flux is completely dominated by the contribution of evanescent 
waves in p polarization II36II37II which depends on the imaginary part of the static reflection 
coefficient which reads for medium i (e,- — 1) / (e + 1). The frequency where £; is approaching - 
1 corresponds to the one where there is the largest density of states of surface polaritons close to 
the material. What appears here is that the contribution of the polaritons of a material is reduced 
when the material is at high temperature. Moreover, one notes that Si02 and SiC do not behave 
in the same way. While the Si02 resonance contribution is weakly attenuated, the one of SiC is 
much more decreased by almost one order of magnitude above 1000 K. This explains why the 
heat flux is larger when the SiC is at low temperature. Moreover rectification does not change 
below a certain distance which correspond to the one where the flux is completely dominated 
by the p contribution of the evanescent waves. It has been shown in several papers 
that this flux behaves as 1/d'^ multiplied by the product of the imaginary parts of the reflection 
coefficients. When the distance decreases, both fluxes increase following the same distance 
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Fig. 4. Computed spectral radiative heat transfer between two planar interfaces made of 
SiC and Si02 in the direct and the reverse situation for 4 different separation distances: 
Si02 is at 297 K and SiC is at 1470 K in the first situation (plain) whereas Si02 is at 1470 
K and SiC is at 297 K in the second situation (dashed). 














scaling law. Therefore, rectification becomes distance independent and saturates. 

5. Conclusion 

We have presented in this paper calculation of the radiative heat flux between two materials 
SiC and Si02 taking into account the real temperature variations of their optical properties. To 
do so, we have used measurements of the optical response of both Si02 and SiC at different 
temperatures from the ambient one up to almost 1500K. We show that when the temperature 
difference between the two samples is large, there is little difference in the far held between 
the fluxes calculated in a situation where one of the materials is at 297 K and the other one 
at higher temperature and another situation where the temperatures are reversed. We see that 
at any temperature difference, rectification increases between 200 nm and 10 nm to reach a 
value which is crescent with the temperature difference. Using the data measured here, one 
has shown that this rectification can reach values up to 0.7, values that are as high as the one 
obtained with phase change materials 02OII39I . These properties could be exploited in thermal 
management. Thermal radiative rectification could for example be adjusted by controlling the 
separation distance between the materials and their temperatures. 
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